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Abstract 

We introduce and study one parameter family of integrable quantum field theo- 
ries. This family has a Lagrangian description in terms of massive Thirring fermions 
ip,tp^ and charged bosons XjX of complex sinh-Gordon model coupled with BCn 
affine Toda theory. Perturbative calculations, analysis of the factorized scatter- 
ing theory and the Bethe ansatz technique are applied to show that under duality 
transformation, which relates weak and strong coupling regimes of the theory the 
fermions transform to bosons and x^X ^^'^ vive versa. The scattering ampli- 

tudes of neutral particles in this theory coincide exactly with S-matrix of particles 
in pure BCn Toda theory, i.e. the contribution of charged bosons and fermions to 
these amplitudes exactly cancel each other. We describe and discuss the symmetry 
responsible for this compensation property. 



^On leave of absence from L. D. Landau Institute for Theoretical Physics, ul. Kosygina 2, 117940 
Moscow, Russia 
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1 Introduction 



Duality plays an important role in the analysis of statistical and quantum field theory 
(QFT) systems. It maps a weak coupling region of one theory to a strong coupling region 
of the other and makes it possible to use perturbative and semiclassical methods for the 
study of dual systems in different regions of the coupling constants. For example, the well- 
known duality between the sine-Gordon and massive Thirring models [|I| plays a crucial 
role for study of many two-dimensional quantum systems. The phenomenon of electric- 
magnetic duality in four-dimensional gauge theories, conjectured in 0| and developped in 
opens the possibility for the non-perturbative analysis of the spectrum and the phase 
structure in supersymmetric Yang-Mills theory. 

Known for many years the phenomenon of duality in QFT still looks rather mysterious 
and needs further study. This analysis essentially simplifies for two-dimensional integrable 
relativistic theories. These QFT's besides the Lagrangian formulation possess also an 
unambiguous definition in terms of factorized scattering theory (FST). The FST, i.e. the 
explicit description of the spectrum of particles and their scattering amplitudes, contains 
all information about QFT. These data permit one to use non-perturbative methods for 
the analysis of integrable QFT and makes it possible in some cases to justify the existence 
of two different (dual) representations for the Lagrangian description of the theory. An 
interesting example of duality in two-dimensional integrable systems is the weak coupling 
- strong coupling flow from the affine Toda theories (ATT) to the same theories with 
the dual affine Lie algebra IQ. The duality in rank r non-simply laced ATTs coupled 
with massive Thirring model was studied in [0]. It was shown there that dual theory can 
be formulated as the non-linear sigma-model with Witten's Euclidian black hole metric 
IP (complex sinh-Gordon theory) coupled with non-simply laced ATTs. Lie algebras of 
these "dual" ATTs belong to the dual series of affine algebras but have the smaller rank 
f = r — 1. 

In this paper in section 2 we consider one-parameter family of integrable QFT, which 
has the Lagrangian formulation in terms of complex fermion field (■?/', 'j/'^), complex bo- 
son field {x,x) and n scalar fields ip = {ipi, ...,{pn). This QFT possesses U{1) ® U{1) 
symmetry generated by fermion and boson charges and Q^. It can be considered as 
BCn ATT coupled with massive Thirring and complex sinh-Gordon (CSG) models. In 
the weak coupling region this QFT admits a perturbative analysis. There the spectrum 
of particles, besides the charged fermions (^,^^) and bosons (x, x), contains the scalar 
neutral particles Ma with the masses characteristic for the BCn ATT. Perturbative cal- 
culations show that classical mass ratios are not destroyed by quantum corrections and 
that charged particles possess non-diagonal scattering. The scattering amplitudes of the 
charged particles can be expressed through the solution of factorization (Yang-Baxter) 
equation. The perturbative analysis together with U{1) ®f/(l) symmetry fix this solution 
up to one parameter, which depends on the coupling constant. 

In section 3 we introduce the external field A coupled with charges and Qy.. We 
use standard Bethe ansatz (BA) technique to find the exact relation between the coupling 
constant and parameter of FST. We show that in the strong coupling regime the behavior 
of fermions ip in the external field is similar to the behavior of weakly coupled bosons x 
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and vice versa. The resulting FST possesses the property of self-duahty together with 
fermion-boson transformation ip ^ x- Another remarkable property of this QFT is the 
exact coincidence of the scattering amplitudes Sab of neutral particles Ma with S-matrix 
of pure BCn ATT. It means that fermion and boson contributions to these amplitudes 
exactly cancel each other. In the last section we describe and discuss the symmetry which 
is responsible for this compensation property. 

2 Integrable deformation of BCn Toda theory and 
factorized scattering theory 

In this section we consider the QFT which can be described by Dirac fermion ip, 
complex scalar field x ^ scalar fields ip = {ipi, ipn) with the action : 
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n-l 



+ -id^^y - 2exp(/?^i) + 2 ^ exp/3iipi+i - ifi) + exp(-2/5y;„) 
^ i=i 



/32 



where g/n = - 4,(1+^2/4^) - 

The last term —^exp{—2(3(pn) in (|1]) plays the role of usual contact counterterm 
which cancels divergencies coming from fermion loop. With this term the action (|Tp has 
the form of BCn affine Toda theory coupled with massive Thirring and CSG |^ models. 
Following the notations of ref. |0] we denote this QFT as BCniipyX: P)- It possesses 
® f/(l) symmetry, generated by the charges 

Q^ = J dx^7o^ ; = / ^^ 2^Xl + (f)^^'^xP) ^ 

The QFT (|T]) is integrable. It possesses the local integrals P, with odd (Lorentz) 
integer spins s. The explicit form of these integrals is not in the scope of this paper. 
Some additional symmetry of i3C„ (■?/', x, /5) model generated conserved charges with half- 
integer spin s = n + I is described in the last section. We also checked at the tree level 
that multiparticles amplitudes factorize into two-particle ones. 

For small f3 we can use the perturbation theory for the analysis of QFT (|l]). Its spec- 
trum contains charged fermions (■?/', ^z^^), charged bosons with mass M and neutral 
particles Ma, a = l,...,n. In one loop approximation the mass ratios in BCniipyX: P) 
theories are not destroyed by the quantum corrections and have the classical values : 

M„ = 2Msin(^); a = l,...,n; (3) 

here and later h = 2n + 1. 

The scattering theory in integrable theory (|ID is completly defined by a two particle 
S-matrix. The non-diagonal scattering is possible only between the particles ipjip^x^X of 
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equal mass. All other amplitudes Sax, ^aip and Sah are the pure phases. The scattering 
matrix of charged particles with U (1) ® U{1) symmetry {ip — >■ e^^ip, x ~^ and C, P, T 
invariance is characterized by the following amplitudes : 



|^(^l)V^(^2) >^n-- 
mi)x{02) >in-- 

mr)ij\e2) >^n 



imx{02mo^) >out + mmo2)x{Oi) >out ; 

= a(^)|x(^2)^^(^i) >out + Pmi^K02)x{0i) >out ; 

= T4e)\^lJ^ie2)^Pie^) >out + R4e)\^{e2)^\e,) > 



out 



\xi0i)xHd2 



T 



+ m\x\02)xiOi) >out +m\x{02)x\0r) ; 

X\02)x{0l) >out + Rxmxi02)x\0,) >out 



> 



out ) 



(4) 



All the amplitudes (||) depend on rapidity difference 9 = 6*1 — 6*2. They satisfy the following 
crossing symmetry condition : 



-- T^{0); R^{m 
7(^); ^^{m-e 



) = R^{6); 



(5) 



The perturbative expansion for these amplitudes in the first order in has the form: 



1 - 



1(3^ 



he-. _ 



Ah 

Rtp{0) = Rx{d) - 



h coth ( 

V 2 ^ 

/icoth (—) + 



sinh 9 
2 



a{e) = 7(^) 



1 - 



2sinh(/i^) 

1(3^ 



2hsmh{e) 



sinh 9 . 

+ 0{(3') 



(6) 



UilTT — 



1(3^ 



4sinh(/i^/2) 



+ 0{(3'). 



Factorization property imposes non-trivial limitations to the scattering amplitudes. 
They should satisfy the functional Yang-Baxter (factorization) relations. There are two 
types of C, P, T invariant solutions of Yang-Baxter equations with f/(l) ® f^(l) symmetry. 
The first one corresponds to the case = and is expressed through the direct product 
of two sine-Gordon S-matrix. It follows from eq. (H) that in the BCn{ip,Xy P) niodel is 
realized another case ^ S^. The solution of the factorization equation with this 
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property possesses one (besides scale of 6) arbitrary parameter x and has the form : 

_ smh{Xe + tx7T) _ sinh(Ag - zxtt) 

^^^'^^ sinh(A0) - sinh(A^) 

T.W^- ^"^'^\^-,7V „W; ^^(,)^ cosh(AW ^^(,)^ 
cosh^Ab') cosh(At^) 
2? sin TTT 

Pie) = 6ie) = -t^^Uoy, m = m=p^^Y^iey / = i. 

smh(At^) cosh(At^) 
The solution is consistent with crossing symmetry (^) only if 

A = ^ = r^+^, P=(-ir^ (8) 

Function Yn{9) satisfies the unitarity and crossing symmetry relations : 

sin ^ TTT 

Yn{e)Y4^e)= ■ ,2,^ ■ 2 ; Y^ie) = Y^i^n -e). (9) 

smh t7 + sm irx 

The minimal solution of eqs. (H) is known |^ and has the form : 

YM - Rmi^^~o), m-\i F,{e,m{OA) ' 

In particular the amplitudes and S^^- which are the pure phases can be represented as : 

duo sm{uj9) sinh(7rxco'//i) cosh[7rc<j(/i + 2(1 — x))/2h] 



oo 



- S^{e) = exp 
-Sxi^) = exp 



oo 00 sin(7rc(j//i) cosh(7rco'/2) 

5° duo sin(cj^) sinh(7r(l — x)uo/2h) cosh[7rci;(/i + 2x)/2h] 
oo uo sin(7rc<j//i) cosh(7rto'/2) 



; (11) 



We note that function Yn{9) is invariant under the transformation x 1 — x. S-matrix 
(|^) is invariant under this transformation together with fermion-boson transformation 

For small x function Yn{9) has the expansion : 
We can see that small x expansion of amplitudes (|^) coincides with perturbative expansion 

m) if 



X 



^{l + OiP')). (13) 



Air 
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The amplitudes (J^) possess poles at the physical strip < Im6 < n which are located 
at the points 9a, where : 

ea = t7r{l-—), a = l,...,n. (14) 

These poles correspond to the neutral bound states Ma with masses (|^). The scattering 
amplitudes Sa^p, Sax, ^ab including these particles can be obtained by the usual fusion 
procedure. For this purpose it is convenient to represent the particles Ma in the form : 



+ (-1)" 



(15) 



and use for the particles ip{6), x(^) the commutation relations with S-matrix In this 
way we obtain the following expression for the amphtudes Sa^ip, Sa^ '■ 

+ 2/3(0 + |)M^-|). 

In particular the amplitudes 5*1^ = 5*1;^ can be written as : 

S (6) S (6) sinh6' — zcos(7r(l — 2x)//i) sinh6' + zcos(7r//i) ^^^^ 
^"^ sinhO + icos{iT{l — 2x)/h) smh.6 — icos{iT/h) 

To describe all two particle amplitudes we introduce the notations Q : 
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si^Mf + W) r.. iz-l)iz + l) 



sinh(l-ff) ' {z-l + 2x){z + l-2xy ' ' 

then 

Sa^ = SaAO) = f[{^ + 2p-a-l}. (19) 
p=i ^ 

The scattering amplitudes of neutral particles Ma have the form : 

la+bl-l 

Sam= n {V\{h-V\- (20) 

|a-b|+l(step2) 

The amplitudes (p!9D, (pO]) are invariant under the transformation x 1 — x. They also 
are in agreement with first order of perturbation theory if x has the form (|13|). To find 
the exact relation between parameter x and coupling constant (3 in action (|l]) we have to 
use non-perturbative approach to the QFT BCn{ip,Xj P)- 
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3 Non-perturbative consideration 



The QFT (|l]) possesses the symmetry U{1) (g) f/(l) generated by the charges and 
dl) and admits the introduction of the external fields and coupled with these 
charges. For simplicity we consider the configurations with only one non-zero field A^ 
or A^, which we denote as A. The hamiltonian 7^^(?i^) in the external field A has an 
additional term equal to —AQjp{—AQ^) : 

= Ho — AQ^ ; = Ho — AQ^, (21) 

where Ho is the Hamiltonian of the QFT (|l]). 

To find the exact relation between x and f3 we calculate the specific ground state 
energy SjIj{A){S^{A)) in the limit A oo from the Hamiltonian and from S-matrix data. 
The calculation of these asymptotics from the Hamiltonian (pT|) follows exactly the lines 
of ref. 0, where similar calculations were done, so we reproduce here the result : 

t^{fj, A oo) = , (22) 

Ztt 

^x(A A^oo) = — . 

We calculate now the same values from the S-matrix using the BA approach (see for 
example refs. 0). We consider the case corresponding to the ground state energy £^ 
taking into account that for function all consideration differs only by the notations 
[ip — > x)- Due to additional term —AQ^ every positively (negatively) charged particle 
'il){0){;ip\0)) acquires the additionnal energy A{—A). For A > M the ground state con- 
taines a sea of positively charged particles '?/'(6'), which fill all possible states inside some 
interval —B < 9 < B. The distribution e^{9) of particles within this interval is deter- 
minated by their scattering amplitude S^{6). The specific ground state energy can be 
expressed through the function e^{6) as : 

S^A) - £^0) = -— / coshee4e)d9, (23) 
Ztt J-b 

where non-negative function e^(^) satisfies, in the interval —B < 9 < B, the BA equation: 

k^{e - e')e^{e')de' = a - Mcosh^, (24) 

J — B 

where the kernel K^{9) in (^) is related to the ipip scattering phase by : 

k^{e) = 5{e)-^.^\ogS^{e), (25) 

2m au 

and the parameter B is determinated by the boundary conditions e^(±i?) = 0. 

The Fourier transform K^{uj) of the kernel (|25| ) can be obtained from eq. (0) and 
has the form : 



sinh 



T!:uj{l — x)/h cosh 7Tuj{h + 2x)/h 



cosh(7ru;/2) sinh{7ruj / h) 



(26) 
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For the function S)^{(3,A) we obtain exactly the same equations with function e^{6) sat- 
isfying eq. (p4D, where kernel K^{0) is related to the amplitude S"^ by the eq. (pS]). The 
amplitudes and (|ll]) are connected by the transformation x — 1 — x. It means 
that: 



KM 



sinh 



TTUJx/h 



cosh 



7rcj(/i + 2(l-x))//i 



cosh(7rct;/2) ?,mh.{T:uj / h) 

The main term of the asymptotics of the function £^ at A ■ 
explicitly through the kernel K^i^uj) {Ky^{uj)) by the relation |Ty] : 

A^ A^ 



(27) 

oo can be expressed 



8^{(3,A 
E^{f5,A 
Comparing eqs. (|2^) and 



oo) 



oo 



27rK^(0) 
A^ 



27r(l 
A^ 



X] 



(28) 



a;(/3) 



27rKx(0) 2'Kx' 
we find the exact value for parameter x(/5) 
, /47rs 



47r + /32 



(29) 



The term in eq. (23) is the bulk vacuum energy of QFT (0). It can also be expressed 
through the kernel K^{ijj){K^{ijj)) by the relation : 



^0 



K^^y^{uj) cosh ■ 



sinh7r//i 



8 sin(7ra;//i) sin(7r(l — x)/h) 



(30) 



The bulk vacuum energy is symmetric under the transformation /3 — 47r//5. It follows 
from the eqs. (P^D, ( P7| ) and ( pUj ) that the transformation relates the function £^ and : 



£^{(3,A)=£^{j,A). 



(31) 



It means that strong coupling behavior (/5 >> 1) of the fermions ip^ip'^ in the external 
field A coincides with weak coupling behavior (/3 << 1) of the bosons x? vice versa. 

In the weak coupling region the behavior of functions £^ and is rather different. In 
this limit K^{6) = 1 + 0(/3^) and function £^{A) can be easily calculated. As function of 
j3 it has the smooth behavior at /3 << 1 and can be written in parametric form as : 



^V,(v4)-^(0) 
2A) 



u{l - 



l-2u + 2m(1 - u) loj 



u 



u 



u) 



(32) 



The kernel K^{6) in the weak coupling limit is not trivial and has the form : 



cosh 



nuj{h + 2)/h 



Ah cosh(7ru;/2) sinh(7ru;//i) 



(33) 
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The BA equations with the kernel ( P5D were studied in [0, where the function Sy^{A) 
corresponding to this kernel was calculated. It can be also written in the parametric 
form: 



\ - 2u - hu(\ - u)[\ - ' ^ 



M\2 ^ , U \2/h 



U 



+ 0{l] 



- 



2AJ J ■ <^^' 

It was shown in [0 that function ( ^4|) coincides with the classical minimum of the 
functional (^1]), where in the weak coupling limit we can neglect the fermionic terms. 
This gives us an additional test relating QFT (|lD and FST (0). The singular behavior 
(0(l//3^)) of function reflects the instability of weakly coupled bosons with respect 
to introduction of external field. The threshold behavior {A = {A — M)/M << 1) of 
functions (^) and p^ ) is also rather different. The first function has there singularity 
~ A^/^ characteristic for fermionic particles. The second one possesses there behavior 
~ characteristic for the weak coupling limit of bosonic theory. It follows from eq. ([31]) 
that properties of particles ip and x change drasticaly under the flow from weak to strong 
coupling (fermions transform to bosons and vice versa). 

The scattering theory of charged particles ip and x (0) is invariant under the transfor- 
mation (3 — > An/ [3 together with fermion-boson transformation ip ^ x- The amplitudes 
Saip and Sax (HD invariant under this transformation. The amplitudes Sab{0) of 
scattering of neutral particles Ma{6) are also self-dual. The remarkable property of the 
amplitudes Sab{0) is that with function defined by eq. ( pQ]) they coincide exactly 
with scattering matrix of particles in pure BCn ATT proposed in ref. IQ. The pure BCn 
ATT can be obtained from the QFT (|I]) by the reduction ip = x = 0- It means that all 
contributions to the amplitudes Sab of neutral particles coming from charged particles x 
and ijj exactly cancel each other. This cancellation (which can be checked in perturbation 
theory) can not be occasionnal. It should be a symmetry responsible for this compensation 
property. We discuss this symmetry in the next section. 



4 Concluding remarks 



The symmetry responsible for the exact compensation of fermion and boson con- 
tributions to the amplitudes of neutral particles should relate ip X particles. It is 
possible only if this symmetry is generated by the conserved charges with half-integer 
spins. It should also be consistent with the FST (0). One can check that this scattering 
matrix commutes with the symmetry algebra T„, generated by the charges (5±(Q±) with 
(Lorentz) spin s equal to h/2 = n + 1/2 (— n — 1/2), and "fermion number" F. The 
charges Q±, Q_^, F possess the following commutation relations : 



Q 



qI 



{Q+,Q_} = {Q- 



F,Q± =±Q± 



± 



= tQ±; 



; 



{Q+,Q.} = P,- {Q+,Q_} = Ph, 



(35) 
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where Ph{Ph) is the right (left) component of the local integral of motion with spin h. 
This local charge acts on asymptotic (in, out) states of charged particles ip{6) and x(^) 
with the eigenvalue A^(6') (A {6)) : 

\{e) = (Mexp^)^/2 . ^ (Mexp 0)^^/1 (36) 

The action of the operators Q±{Q±) and e^'^^ on the one-particle states are defined as : 

Q^\m >= mmo) > q+mo) >= homho) > 

Q.\^^{0)>=x{e)\x{e)>; Q^\m>=m\m>; (37) 

e'"^lx(^) >= e^"1x(^) > ; e*"^|V^(0) >= -e-*"^|V^(^) > ; 
e'^^lW) >= e"*"1x(^) > ; e'^'^l^^ie) >= -e*"^|^^(^) > . 

The operators act at the same way as Q± with the substitution A — > A. The action of 
the charges Q±{Qj.) on many particles sates can be defined from the following co-product 
rules : 

A(Q±) = Q±®/ + e±'"^®Q±, 

=^T®^ + e^"''®^T- (38) 

Using eqs. (|37D, ([38| ) and commutativity condition of algebra T„ with S-matrix (^, we 
can derive all the ratios of the amplitudes (^. 

The eigenvalues of the operator exp(i7rF) on the states \x > >) move from 1 (—1) 
at /3 = 0, to —1 (1) at /? = oo. Particles x ('^) have fractional values of fermion number 
equal to x{f3) {l—x{(3)). At the "self-dual" point = An the fermion numbers of particles 
ip and X s^re equal to 1/2. At this point we have the symmetry ip ^ x- The amplitudes 

and coincide and S-matrix (^ of the QFT BCni'ip, X) V^n) can be expressed as the 
direct product of two-particle S-matrices of sine-Gordon model : 

, , -JV sinh ^ — i sinfTra/Zi) , ^ , ^ , ^ 

^J-Q smh + z sm{na/h) 

where Sn{6) is S-matrix of sine-Gordon model corresponding to the coupling constant 
^sg/^'^ ~ '^^^ FST (39) coincides with S-matrix of Cl^li{ipi^2, VAtt) model consid- 

ered in 0. It means that corresponding QFTs also coincide at the point = An. 

To construct the currents ^^+3/2(^1+3/2) generating conserved charges (5±(Q+) it is 
convenient to rewrite the action (|I|) in terms of other fields. We introduce the scalar field 
$ related with the fields ipjip"^ by the usual bosonisation rules [0, and the fields ^o^^Po 
which describe the dual representation for CSG-model [0, |ll|. In terms of these fields the 
action (0) can be represented as : 

An = J d^x [^{d^%Y + lid^^oY - Mo cos(a'$o)e^'^» (40) 
+ -id,^r - Mo cos(«$)e-^'^" + -id.^f - -JI [ ^ exp /3(^,+i - 
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where parameters a, a', j3' are defined by the relations : 

- = 47r , a'^ - (3'^ = Air , = Air/ (3. (41) 

The first three terms in (^) which we denote as Af^{^Q, (fo) correspond to the first (sigma- 
model) term in the action (|l]). With action Ao-($o,V^o) the fields X:X and XX can be 
represented as : 

/2z7r 27r \ _ / 2z7r 27r n 

X ~ exp (^—$0 - —vo) ; X ~ exp - —$0 - — </^oj ; (42) 

XX - const ~ exp - ^V^o) = exp(-/5v?o)- 

The action (^) is not available for the ordinary perturbation theory in coupling constant 
P, however with this action QFT BCni^p, X? VAtt) can be treated as perturbed conformal 
field theory. Using this approach we introduce the fields (/>o(0o), 0(0) which are the right 
(left) chiral components of fields $0, ^, and right (left) derivatives d = do + di {d = 
do — di). The straightforward calculation shows that spin n + 3/2 (— (n + 3/2)) currents 

^n+3/2(-^i+3/2); which generate conserved charges Q±{Q^)can be written in the form : 

^n+3/2 = exp ± —(pjiKd - dipn)...{Kd - dipi){Kd - difo) exp ± -^0oj, (43) 

where k = (-^ + — 

^An j3 

The currents can be obtained from the eq. (|^) by the substitution 



n.+3/2 can ue ouiaineu iiom me eq. uy ine suuslilulioii cpo ~^ •f'O) 

— s> 0, 9 — >• 5. 

The fermion number F can be expressed through the charges and (^ as : 

F = x(/3)Q^ + (l-x(/5))Q^ . (44) 

Together with charges Q± and it generates the algebra T„. It follows from the rep- 
resentation (|15D for the neutral particles M^iff) and eq. (|36D that all these particles are 
annihilated by the local conserved charge '■ 

P^\MaAOl)...M,Adn) >rn(out)= 0- (45) 

The neutral boson sector of the theory which contains only asymptotic particles Ma{9), 
is defined by the kernel of the conserved charge Ph- The restriction of the scattering theory 
to the neutral sector DT, defines self-consistent FST. The QFT corresponding to this FST 
possesses the Lagrangian description with the local action, which can be obtained from 
the action (P by the reduction ip = x = ^■ 

At the end we note that besides this reduction, the QFT (0) possesses several other 
integrable reductions. For example : x = 0, '0 7^ 0; -0 = 0, x 7^ 0; 7^ 0, x is real; 
X 7^ 0, 0^ is Majorana fermion; x is real, is Majorana fermion and so on. The FSTs 
corresponding to these reductions are described in refs. ||^, 0, |]. 
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